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Abstract 

By means of micromagnetic spin dynamics calculations, a quantitative calculation is carried out 
to explore the mechanism of exchange bias (EB) in ferromagnetic (FM) /compensated antiferro- 
magnetic (AFM) bilayers. The antiferromagnets with low and high Neel temperatures have been 
both considered, and the crossover from negative to positive EB is found only in the case with low 
Neel temperature. We propose that the mechanism of EB in FM/compensated AFM bilayers is 
due to the symmetry broken of AFM that yields some net ferromagnetic components. 
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Exchange anisotropy was first discovered in 1956 by Meiklejohn and Bean[lj], who found 
that the hysteresis loop of Co/CoO after cooling in a magnetic field was no longer centered 
at zero field (H=0) but was shifted along the field axis. The shifted direction was found to be 
opposite to the applied magnetic field (negative exchange bias Hb < 0) and the magnitude 
of this shift is known as exchange bias (EB). It was subsequently established that this 
might be a general phenomenon for any ferromagnet (FM) / antiferromagnet (AFM) systems 
cooling in an applied magnetic field (cooling field) from above the Neel temperature(TTv) 
of the AFM, with the FM Curie temperature(T(7) greater than TV- In recent years, since 
the phenomenon of exchange bias has become the basis for an important application in 
information storage technology 0, tremendous efforts have been made for exploring the 
mechanism ||, Q. 

Meiklejohn and Bean originally suggested that exchange bias was a consequence of the 
presence of interfacial uncompensated AFM spins. In view of this argument, a natural 
question to ask is whether the exchange bias also exists in a FM/ compensated AFM system. 
Surprisingly, in a compensated Fe/FeF2 bilayer system, Nogues et al. observed not only 
the usual negative exchange bias but also an unexpected positive exchange bias (Hb > 0) 
under large cooling fields p. 

Several important theories have existed to study the exchange bias in compensated 
AFM. Koon|Bj presented a microscopic explanation of EB due to a irreversible AFM do- 
main wall, and found a perpendicular orientation between the FM/AFM axis directions, 
namely spin-flop state. With consideration of magnetostatic interactions in this spin-flop 
state, Schulthess et al. |7j obtained the opposite results, i.e., not EB but a large uniaxial 
anisotropy, and attributed EB to the interfacial defects. Unfortunately they did not further 
show the EB-magnetostatic interactions phase diagram, in other words, EB should gradu- 
ally change with magnetostatic interactions. Hong|| argued that interface spin configuration 
persisted after cooling below Tjv, and negative/positive bias respectively corresponded to 
parallel/perpendicular easy axes of FM and AFM. Kiwi et al.@ suggested a canted AFM 
spin configuration frozen into a metastable state, and proposed the incomplete FM domain 
wall model to explain positive exchange bias. However, the former two theories were carried 
out with micromagnetic calculations without consideration of the cooling field; the later two 
theories pointed out the cooling field without micromagnetic calculations, and were lack 
of much more detailed and sufficient microscopic information. Up to now, exchange bias 
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mechanism is still controversial. 

In this paper, based on the assumptions that an antiferromagnetic interface coupling be- 
tween FM/AFM is responsible for exchange bias in FM/compensated AFM ||, |TIJ- |I2| , and 



that the biased hysteresis loop is basically determined by the spin configurations in the under- 
lying antiferromagnetic layer after cooling [TB|]-[T5|], we carry out micromagnetic calculations 



using spin dynamics to explain the mechanism of EB in the FM/compensated AFM systems. 
Physically the key point different from previous Koon's micromagnetic calculations 0, ad- 
dresses the cooling field during cooling process. We succeed to reproduce both the negative 
and positive EB effects. Qualitatively speaking, it is a competition among (i) the cooling 
magnetic field (ii) the interface coupling of FM/AFM, and (iii) the spin-spin interaction and 
anisotropy of AFM, that eventually determines the spin configurations in AFM during the 
cooling process. For an AFM with weak spin-spin interaction (low TV), the spin configura- 
tion of AFM at low cooling magnetic field is dominated by the AF type interface coupling 
of FM/AFM. Therefore, the initially compensated AFM layers especially the interface AFM 
layer becomes weekly uncompensated, resulting in a net ferromagnetic component opposite 
to the cooling field (or the magnetization in FM). The hysteresis loop is then measured at 
low temperature after removing the cooling field, while the spin configuration in AFM is 
frozen. Similar to the arguments given by Meiklejohn and Bean, it can easily be deduced 
that the broken symmetry of AFM in this case favors the negative exchange bias. However 
on the other hand, if the cooling field Hcf is large and becomes to dominate, then a net 
ferromagnetic component along the cooling field is expected. Because of the AFM type 
interface coupling of FM/AFM, it turns out that the broken symmetry of AFM in this case 
favors the positive exchange bias. On other hand, for AFM with strong spin-spin interaction 
( high Ttv ) only negative EB can be found in a reasonable high Hcf- The quantitative 
results are given in the following to reveal in details how these different terms affect the 
broken symmetry of AFM layers and its correlation with the exchange bias. 
Our model Hamiltonian is 

H = Ua-a + Hp__p + H A _ F , (1) 
where Ha~a is the part of AFM layers, Hp_p and Ha~f the FM and the interface coupling 
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between AFM and FM layers. They are 

Ha-a = E J A-ASi • Sj - D A J2 (S*) 2 - HcfY, (9au b ) S?, (2) 

<i,j> i i 

Hf-f = — X/ JF-F$i • Sj — Hcf2_j (9fub) Sf, (3) 

<i,j> i 

Ha-f = Ja-fSo, ■ S/, (4) 

«*,/> 

where g^, g^, «b, Da and i/o.F denote AFM Lande factor, FM Lande factor, Bohr magneton, 
antiferromagnetic anisotropy and cooling field in parallel with AFM anisotropy, respectively. 
The exchange coupling among spins is considered for nearest neighbor sites only. The 
subscripts a and / are associated with AFM and FM respectively. It is noticed that the 
anisotropy of FM layer is neglected based on the fact that most experiments did use the soft 
ferromagnets. The dipole-dipole interactions in the system are not considered here, since 
it affects only quantitatively rather than qualitatively on the symmetry broken of AFM. 
As the previous models, we also assume that J a-f ~ J a- a @, 0- The Neel temperature 
increases monotonically with J a- a, thus the interface coupling is stronger in the FM/AFM 
systems with higher Tjv and vice versa. 

Now we calculate the EB by the following spin dynamics approach P], i.e., the local ef- 
fective field is determined from the gradient of the energy, H^ ff = — ^ s , and {Sj} is 
required to satisfy the Laudau-Lifshitz equation of motion with the Gilbert-Kelley form for 
the damping term: J^-Sj = gu B Si x (Hf ff — r/^S^, where rj denotes the damping param- 
eter. This damping term is phenomenological and is included to remove the energy from 
the system and to ensure that the magnetic system is in a stable or metastable equilibrium 
after sufficient iterating calculating steps. A lattice with 50 x 50 x 2(FM) and 10 layers of 
AFM is used in our calculation. In the beginning, the temperature T of the system is set 
at T c > T > Tjv- Therefore, the initial spin configuration in our spin dynamics calculation 
is such that the spins are randomly arranged for AFM but are ferromagnetic arranged for 
FM. A cooling field is then applied along the x direction that is also the easy axis of AFM. 
Meanwhile it is known that the spins of FM will be easily aligned according to the cooling 
field. After taking a long step of spin dynamics calculation, a stable state of FM/AFM 
under the cooling field is finally approached. Then the system is cooled down to the low 
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temperature, we switch off the cooling field and start to do simulation of hysteresis loop of 
FM layers while the spin configuration of AFM is fixed. 

Fig.l shows EB as a function of cooling field for a AFM/FM system with low TV such 
as FeF 2 (TV ~ 78 AK) or MnF 2 (2V ~ 67.3f^). In doing this we set the parameters in 
Hamiltonian as qa = 9f = 2.0, Jf-f = lOmev, J a-a = 0.8mev, J a-f = Ja-a/2 = OAmev 
and Da = OAmev per site0, |J, J7|. As a natural output from the calculation, it is indeed 
observed in this figure that the exchange bias Hb changes sign from negative to positive as 
the cooling field increases, and a crossover field H cross is found at about 3.7T. Dashed and 
solid lines in the inset show the negative and positive loops at 2KOe and 7T respectively. 

For systems with low TV, i.e., weak spin-spin interaction Ja-a in AFM, the spin config- 
uration of AFM at low cooling magnetic field is dominated by AFM type interface coupling 
of FM/AFM. In this case, the symmetry broken of compensated AFM layers appears. Some 
net ferromagnetic component along the — x axis is expected, which means that the broken 
symmetry of AFM in this case favors the negative exchange bias. However for higher cool- 
ing magnetic field Hqf, the cooling field becomes to dominate the broken symmetry so that 
a net ferromagnetic component along the positive x direction is expected, i.e., the broken 
symmetry of AFM in this case favors the positive exchange bias. 

For a quantitative description of the FM components in AFM, We define the ferromag- 
netic component in the n-th AFM layer 

S x (n)=J2SL(i)/N n , (5) 

i 

where S x a (i) is the AFM spin at site i of the n-th layer, N n is the number of lattices in the 
n-th AFM layer, while the first layer is defined as the interface layer of AFM. This quantity 
S x (n) describes the degree of symmetry broken in each layer of AFM. It is found from our 
calculation that the ferromagnetic components are layer dependent. As expected, it will be 
larger when the layer is near the interface, and become smaller when the layer is far from the 
interface. In Fig. 2 the ferromagnetic component of the interface AFM layer S x (l) is shown 
as a function of Hqf, using the same parameters as obtaining Fig.l. Similar results as 
Fig.l are found, that S x (l) is negative at the beginning when the cooling field is small, then 
reaches to zero at a critical field (~ 37KOe) and finally becomes positive as the field further 
increasing. Since the first AFM layer should be the most important one in the interface 
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coupling, it is reasonable to see that the ferromagnetic component of the interface AFM 
layer should be responsible for the EB effect. The inset of Fig. 2 gives the layer dependent 
ferromagnetic components when the cooling field is fixed at 2KOe. The oscillation is caused 
by the antiferromagnetic exchange interaction between the layers of AFM. 

Fig. 3 shows the same relationship but for a AFM/FM system with higher T/v such as 
FeMn (~ 500K). The parameters used here are Jf-f = lOmev, J a- a = 5mev, J a-f = 
J a- a = 5mev and Da = 3mev per site |§ . In this case, it is found that Hb is always negative 
and changes a little when Hcf is ranged from 2T to 7T. This result also agrees qualitatively 
with experiment 0. Inset (a) of Fig. 3 shows the magnetization loop at 3T cooling field. In 
fact that J a- a is large when the AFM layer of the system has high T/v, thus the interface 
coupling Ja-f also becomes large (Ja-f ~ J a- a), then the AFM type interface coupling 
controls the symmetry broken of AFM. In this case, if Hcf is reasonable high (2T — 7T) 
but not too high, Hb is found to be always negative. This can explain that the positive EB 
was reported only in the FM/AFM thin films with low T/v@, [10] - [12 . 



One distinguished feature of Fig. 3 is that: a tip EB is found at cooling field H t i V indicated 
by arrow, and this tip also is clearly shown in Inset (b). As previous mention, the preceding 
discussions are subjected to both cooling field and applied magnetic field parallel with AFM 
easy axis. In fact, the orthogonal FM/AFM spin configuration similar to Koon's conclusion 
can also be recovered with zero or smaller cooling field for stable spin configuration. With 
the increasing cooling field from zero to Hu p , the FM spins will gradually rotate direction 
from the perpendicular to parallel to the AFM easy axis during cooling process, and the 
interface coupling contribution to the negative bias will enlarge and nearly saturate at Hu p . 
On the other hand, with the increasing cooling field above H tip , the cooling field contribution 
to the potentially positive bias will raise, in other words, contribution to the negative bias 
will lessen. Thus for cooling field parallel with AFM easy axis situation, there exists a tip 
EB associated with cooling field Hu P - 

In summery, micromagnetic spin dynamics calculations are carried out to explain the 
mechanism of EB in the FM/compensated AFM system. Different from the previous micro- 
magnetic calculations, we address the key role of cooling field. Some important experiment 
results, such as the cooling field dependent transition from negative to positive EB in the 
AFM/FM layers with low T/v, can be reproduced. It is proposed that the symmetry broken 
of AFM plays a key role in explaining the exchange bias. 
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FIG. 1: Exchange bias Hb as a function of cooling magnetic field Uqf for FM/AFM with lower 
Tj\r {FeF-i). Dashed and solid lines in the inset show the negative and positive magnetic loops at 
2KOe and 7T cooling field, respectively. 
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FIG. 2: The ferromagnetic component along the x axis as a function of cooling field for the interface 
AFM layer. Inset shows the layer dependent ferromagnetic components of AFM layers at 2KOe. 
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FIG. 3: The relationship between exchange bias Hb and cooling field Hcf for FM/AFM with 
higher Tjv (FeMn). Inset (a) presents a magnetic loop at 3T cooling field. Inset (b) shows clearly 
the EB around cooling field H t i p indicated by arrow in the figure. 
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